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摘 要 :根系 是 评价 植被 生态 服务 功能 的 关键 基础 要 素 ,但 快速 精确、 无 损 确 定 根 系 的 测量 技术 和 方法 是 目前 生 


态 系统 评 佑 中 的 瓶颈 。 


探 地 雷达 是 一 种 高 效 无 损 的 地 球 物理 学 技术 ,可 以 在 无 损 状态 下 获取 土壤 中 的 根系 信息 。 


但 是 探 地 雷达 检测 和 识别 根系 的 精度 受 土壤 含水 量 、 根 系 含水 量 、 根 径 大 小 、 根 系 埋藏 深度 等 诸多 因素 影响 ,导致 
其 在 野外 根系 探测 中 适用 性 受 限 ,为 了 探究 土壤 含水 量 对 探 地 雷达 探测 根系 精度 的 影响 ,本 人 研究 采用 野外 预 埋 根 


系 的 控制 实验 ,根据 探 地 雷达 波 速 .振幅 和 根 点 反射 系数 的 变化 ,分 析 了 不 同 土壤 含水 量 条 件 下 根 点 识别 率 及 根 点 
距离 均 方 根 误差 。 结 果 表 明 :(1) 根系 探测 中 , 探 地 雷达 波 速 和 振幅 是 判断 土壤 含水 量变 化 的 重要 参数 ;(2) 随 着 
土壤 含水 量 的 增 大 , 探 地 雷达 波 速 减 小 ,雷达 振幅 趋 于 平缓 ;(3) 不 同 土 层 深度 上 随 着 根系 直径 的 增 大 , 探 地 雷达 波 


速 增 大 ,雷达 振幅 趋 于 激烈 ;(4) 根 点 识别 率 与 土壤 含水 量 成 负 相关 (P<0.05) ,土壤 合 水 量 为 15%~25% 时 探 地 雷达 
对 活 根 的 识别 效果 最 佳 。 本 研究 表明 探 地 雷达 可 以 作为 植物 根系 生物 量 无 损 和 快速 测定 、 评 估 的 方法 ,但 在 利用 
探 地 雷达 测定 土壤 中 根系 时 ,应 在 土壤 含水 量 相对 较 低 的 时 间 段 进行 。 

关键 词 : 探 地 雷达 ， 根系; 土壤 含水 量 ， 无 损 探测 ， 亚 西 黄土 区 


探 地 雷达 (Ground Penetrating Radar, GPR ) 通 过 
天 线 将 电磁 波 发 射 到 地 下 ,利用 不 同 介质 的 电磁 特 
性 差异 ,产生 不 连续 的 反射 和 散射 ,从 而 定性 和 定 
量 地 辨别 地 下 目标 物 ,实现 对 目标 物 的 无 损 、 快 速 
和 精确 探测 ,因此 , 探 地 雷达 作为 一 种 高 效 无 损 的 
地 球 物理 学 技术 ,被 广泛 应 用 于 物探 领域 的 探测 、 
成 像 和 特征 识别 等 物探 领域 中 ”。 

根系 在 植物 生长 发 育 过 程 中 从 土壤 中 汲取 水 
和 营养 物质 ,并 起 着 支撑 和 固定 植物 体 的 作用 ,更 
是 防治 水 土 流失 ,保育 土壤 的 重要 带 官 。 植 物 的 各 
种 生态 服务 功能 均 与 根系 生物 量 密切 相关 ,因此 ， 
根系 生物 量 的 准确 测定 是 科学 评 佑 条 林 生态 服务 
功能 的 基础 和 关键 。 但 由 于 土壤 不 透明 属性 的 限 
制 ,埋藏 于 土 层 中 错综复杂 的 植物 根系 难以 直接 观 


收 稿 日 期 : 2023-09-11; 修订 日 期 : 2023-11-30 
基金 项 目 : 国家 重点 研发 计划 (2022YFE0104700) 
作者 简介 : 
通讯 作者 : 张建军 . E-mail: zhangjianjun@bjfu.edu.cn 


察 与 测量 ,人 工 调查 植物 根系 耗 时 费力 ,调查 结 
还 具有 随机 性 。 正 是 由 于 根系 调查 方法 具有 很 大 
的 局 限 性 和 挑战 性 ,限制 了 根系 研究 的 发 展 ,进而 
影响 了 对 森林 生态 系统 生态 服务 功能 的 科学 评价 一 。 
已 有 研究 表明 相 邻 介质 之 间 的 介 电 常 数 差异 ， 
决定 了 探 地 雷达 电磁 波 反射 能 量 的 强度 ” ,植物 根 
系 含 水量 高 于 周围 土壤 含水 量 ,根系 与 土壤 之 间 的 
介 电 党 数 人 存在 明显 差异 ,这 为 利用 探 地 志 达 探测 根 
系 提供 了 可 能 性 ””。 近 年 来 探 地 雷达 已 被 应 用 于 
植物 根系 的 无 损 性 探测 ,与 传统 探测 方法 (如 挖 根 
法 AGATA 土壤 剖面 法 ) 相 比 , 探 地 雷达 在 探测 时 
不 损伤 植物 、 不 扰动 土壤 ,在 较 短 的 时 间 内 能 够 完 
成 对 野外 较 大 范围 的 根系 探测 ,具有 连续 性 和 可 重 
复 性 ,与 微 根 窗 管 法 .核磁 共振 成 像 .高 分 辨 率 X 射 
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线 断 层 扫描 (CT) 和 电阻 率 成 像 技 术 (ERT) 等 根系 
探测 技术 相 比 , 探 地 雷达 有 具有 投入 少 灵活 性 强 , 可 
在 野外 较 大 范围 检测 植物 根系 的 特点 ""。 

Xiao 等 中 利用 探 地 雷达 对 毛竹 (Phyllostachys 
Pubescens) 根 系 进行 研究 后 指出 , 探 地 雷达 是 根系 探 
测 中 的 可 行 方法 。 但 探 地 雷达 对 根系 的 探测 效果 
受 诸 多 干扰 因素 的 影响 ,其 探测 能 力 在 60 cm 土 层 
Py Bean , 且 随 着 土 层 深度 的 增加 而 减弱 。 探 地 雷 
达 对 土 层 中 根 径 较 粗 的 根 识 别 率 更 高 ,可 识别 的 极 
限 根 径 为 5 mm’), WEIS Pe HA! aR Hh FE IA AE TE TE 
定 环 境 下 有 效 识 别 植物 根系 。Bi 等 "和 Sun SF 
究 发 现 基于 深度 学 习 方 法 下 , 探 地 雷达 对 根系 定位 
和 根 径 估算 的 精度 会 更 加 准确 。Dannoura 等 ” 认 
为 根系 与 土壤 含水 量 的 差异 是 影响 探 地 雷达 对 根 
系 检 测 能 力 的 关键 因素 。Bain °° WY WSS AB 4 
土壤 含水 量 与 根系 含水 量 之 间 无 显著 差异 时 , 探 地 
雷达 对 根系 的 识别 效率 显著 下 降 。 但 是 ,对 利用 探 
地 雷达 调查 根系 生物 量 时 最 为 适宜 的 土壤 含水 量 
方面 的 研究 ,以 及 利用 土壤 含水 量 对 探 地 雷达 调查 
结果 进行 修正 等 方面 的 研究 有 待 进一步 深入 。 为 
此 本 人 研究 在 山西 吉 县 森林 生态 系统 国家 野外 科学 
观测 研究 站 的 所 在 地 蔡 家 川 流 域 ,采用 时 外 预 埋 植 
物 根系 的 控制 实验 ,通过 雷达 波谱 图 像 分 析 , 定 量 
人 研究 不 同 土壤 含水 量 对 探 地 雷达 波 速 .振幅 和 根 点 
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识别 率 的 影响 ,以 期 为 提高 探 地 雷达 野外 根系 探测 
精度 和 正确 佑 算 根系 生物 量 提供 参考 依据 A 
黄土 区 典型 人 工 林 水 土 保持 效益 评价 和 碳 汇 评价 
提供 科学 依据 。 


1 实验 与 方法 


1.1 实验 场地 

实验 场地 位 于 山西 吉 县 木林 生态 系统 国家 野外 
科学 观测 研究 站 所 在 地 蔡 家 川 流 域 (图 Ta) ,地 理 坐 
标 为 36°14'27”~36°18'23”N, 110°39'45"~110°47'45"E, 
海拔 为 904~1592 m, 流 域 呈 东 西 走向 ,地 势 西 高 东 
低 ,流域 面积 为 39.33 km’, 流域 长 度 为 12.15 km, JE 
于 黄土 梁 状 丘陵 沟 蜜 区 ,黄土 母 质 ,流域 内 土壤 为 
褐 土 , 土 层 厚 度数 十 米 ,土质 较 均匀 ””。 该 研究 区 
属于 暧 温带 大 陆 性 气候 ,年 平均 气温 10 CC。 平均 降 
水 量 575.9 mm, 多 集中 在 7 一 10 月 ,年 平均 蒸发 量 
1687 mm, 年 平均 日 照 时 间 2563.8 h ,无 霜 期 172 do 
主要 土地 利用 类 型 有 林地 `. 农 地 no SO. BR p A5, 
主要 树种 为 刺槐 、 油 松 、 侧 柏 等 。 
1.2 研究 方法 

实验 材料 油 松 (Pinus tabulaeformis) ,根系 发 达 ， 
适应 性 强 , 耐 干旱 贫 次 ,具有 良好 的 保持 水 土 涵养 
水 源 和 改良 土壤 作用 ” ,是 黄土 高 原 主要 造林 树 
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图 1 样 地 示意 图 
Fig. 1 Chematic diagram of the sample plot 
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种 ,因此 实验 选取 油 松 (Pinus tabulaeformis ) 为 研究 
对 象 。 本 研究 挖 取 不 同 粗 度 的 油 松 根系 作为 实验 
材料 ,将 挖 取 的 油 松 根 截取 成 30 em 的 根 段 作为 样 
根 , 并 在 样 根 两 端 涂 抹 乳 胶 防止 水 分 散失 ,用 游标 
卡尺 量 取样 根 直 径 ,用 电子 天 平 测量 样 根 湿 重 。 
实验 坑 制作 如 图 1b, 在 1990 年 营造 的 人 工 油 松 
林内 挖 1.0 mx0.6 mx0.4 m 的 样 坑 , 将 样 坑 左 边界 和 
上 边界 的 交点 确定 为 坐标 原点 ,原点 垂直 回 土 层 方 
癌 为 Z 轴 ,左边 界 、 上 边界 分 别 定义 为 Y 轴 和 半 轴 。 
ETE HH AY EE AIR EA BER a if Jar VE A 
JR SH. A HE ST EY ESS, S 
400 MHz < FS 800 MHz 的 探 地 雷达 天 线 能 探测 根 
径 大 于 20 mm ,深度 在 2.5 m 以 内 的 根系 。 
样 根 布设 如 图 2 所 示 , 本 研究 将 根 径 2.0 cm 和 
2.5 cm 的 样 根 每 间隔 10 em 放置 于 实验 样 坑内 距 地 
面 10 cem、20 cm,30 cm 土屋 处 ,共计 16 根 ;将 根 径 
4.5 cm 的 样 根 每 间隔 30 cm 放置 于 实验 样 坑 内 距 地 
面 10 cem、20 cm、30 cm 土屋 处 ,共计 3 根 ,所 有 样 根 
放置 均 与 样 坑 上 下 边界 平行 ,分 别 记录 样 根 在 实验 
样 坑内 的 空间 坐标 (X、Y.Z) 后 ,用 过 筛 后 无 其 他 根 
系 的 均 质 土 填 埋 实验 样 坑 。 填 满 后 在 实验 坑 表面 
每 隔 0.2 mm 设置 1 条 长 1m 的 红色 尼龙 绳 作为 测 线 
(图 1lec) ,共计 3 条 , 测 线 设 置 与 样 根 预 埋 方 向 垂直 。 
实验 坑 制 作 好 后 进行 次 水 至 土壤 水 分 饱和 , RF 
实验 坑 地 表 无 积 水 后 ,用 探 地 雷达 (CAS.S500， 
IECAS, China) 主 机 和 频率 为 500 MHz 的 天 线 , 沿 3 
条 测 线 进行 扫描 ,扫描 频率 为 51.2 Hz, 测 量 方式 为 
测量 轮 连续 点 测 。 其 后 每 隔 3 d 对 实验 坑 利 用 探 地 
雷达 进行 一 次 扫描 ,共计 扫描 5 次 ,每 次 扫描 后 使 用 


(a) 土 层 深度 10 cm (b) 土 层 深度 20 cm 
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铝 盒 按照 0~10 cm, 10-20 cm、20~30 cm 分 层 取样 ， 
用 人 烘 干 法 测 得 土壤 含水 量 。 

探 地 雷达 的 探测 原理 (图 3) , 探 地 雷达 通过 机 
载 天 线 将 高 频 电磁 波 从 地 表 以 宽频 短 脉冲 的 形式 
发 射 进 入 地 下 ,经 地 下 界面 和 局 部 目标 体 反 射 后 返 

回信 号 被 天 线 接收 并 记录 ”。 电 磁 波 在 地 下 介质 
中 传播 的 路 径 、 电 磁场 强度 与 波形 与 介 质 的 电 性 、 
几何 形态 等 因素 密切 相关 ”。 因 此 可 以 根据 回 波 
言 号 的 时 延 、 普 特 性 等 参数 ,解释 出 根系 粗 
FE 分布 深度 .介质 结构 性质 及 空间 分 布 等 特征 2 。 
1.3 数据 处 理 

数据 使 用 RadarViewV18.11 和 ReflexW 9. 5 软 
件 进 行 处 理 。 在 RadarViewV18.11 中 进行 背景 直流 
去 除 和 自动 增益 ,然后 进行 F-K 偏 移 和 滤波 处 理 ,最 
后 根据 数据 需求 进行 分 段 增益 ;在 ReflexW 9.5 中 抽 
取 平 均 iH (subtract- mean), Wi IE (static correction ) , 
手动 增益 (y) (manual gain (y) ) , 背景 去 除 (back 
ground removal) , 带 通 滤 波 (band passbutter worth) 。 
在 ReflexW 9. 5 中 得 到 预 处 理 的 雷达 图 像 后 ,进入 模 
块 2D-data annalysis 中 ,使 用 pick the reflector-phase 
follower 进行 相位 追踪 ,保存 数据 后 再 使 用 Create 
LayerShow 功能 提取 雷达 波 速 数据 信息 。 

数据 统计 和 分 析 使 用 Excel 和 Origin 2022 ,图 像 
后 期 处 理 采 用 Photoshop 2022 ,在 SPSS Statistics 24 
中 使 用 方差 分 析 对 雷达 波 速 进行 显著 性 检验 。 

1.4 精度 验证 

本 文 的 根 点 识别 率 参考 张 璐 云 等 ”的 根 点 识 
别 定义 , 即 对 于 特定 条 件 下 所 有 剖面 假设 实测 根 点 
数量 为 内, 某 一 频率 对 应 剖面 的 雷达 图 像 中 提取 出 


(c) 土 层 深度 30 cm 
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注 :实验 区 范围 1.2 mx0.8 mx0.4m ,扫描 区 范围 1.0 mx0.6 mx0.4 m, 采 样 区 范围 1.0 mx0.1 mx0.4 mo 
图 2 样 根 布置 示意 图 


Fig. 2 Sample root layout diagram 
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式 中 :NN 为 匹配 到 的 根 点 对 数 ;% 表 示 第 i 个 实际 预 
埋 根 点 的 横 坐 标 ;表示 第 j 个 雷达 探测 根 点 的 横 坐 
标 ;% 表 示 第 i 个 实际 预 埋 根 点 的 纵 坐 标 ;表示 第 j 
个 雷达 探测 根 点 的 纵 坐标 。 


2 结果 与 分 析 


2.1 土壤 含水 量 对 探 地 雷达 探测 波 速 的 影响 

实验 测定 过 程 中 第 1 d、 第 3 d. 第 6 d、 第 9d 和 
第 12 d 探 测 时 ,0~30 cm 土 层 的 平均 土壤 体积 含水 
量 分 别 为 37.3% 29.2% ,25.396 22.5% Fil 15.6% ,不 
同 含水 量 条 件 下 探 地 雷达 探测 的 平均 波 速 如 表 1 和 
图 4 所 示 。 

不 同 含水 量 条 件 下 探 地 雷达 的 波 速 有 较 明 显 的 
差别 (P<0.05)。 第 1d 探 测 时 ,10 cm,20 cm 和 30 cm 
土 层 处 探 地 雷达 平均 波 速 分 别 为 0.036 mns'、 
0.058 m*ns ! 410.081 m*ns 5*8 12 d 探 测 时 , 10 cm、 


扫描 位 置 /m 


双 程 时 间 /ns 


t 


注 : 图 中 Tx 为 雷达 的 发 射 器 ,Rx 为 雷达 的 接收 器 ， 
1 为 雷达 信号 的 双 程 时 间 。 
图 3 雷达 双 曲 线形 成 示意 图 


Fig. 3 Schematic diagram of radar hyperbola formation 
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可 以 用 公式 (1) 表 示 : 
E 20 cm 和 30 cm + Jz Ab Tg Hl E 3 E HY Ui BE 43 3 OI 
D= zp X 100% (1) — 0054 m- ns” .0.080 m-ns $10.089 m.ns: ,可见 , 随 


计算 雷达 探测 到 的 根 点 和 实际 预 埋 根 点 的 距 
离 均 方 根 误差 (RMSE) ,以 此 作为 探 地 雷达 探测 根 
系 的 精度 参考 ” ,用 公式 (2) 表 示 : 


VY To 


RMSE = N (2) 


着 土壤 含水 量 的 降低 , 探 地 雷达 波 速 逐渐 增 大 , 且 

随 土 层 深度 的 增加 ,雷达 波 速 呈 上 升 趋势 。 

2.2 土壤 含水 量 对 探 地 雷达 探测 波谱 振幅 的 影响 
探 地 雷达 归 一 化 振幅 变化 如 图 5a 所 示 ,5 次 扫 

描 雷 达 振 幅 主 要 集中 在 262.920 Hz~788.762 Hz。 从 

振幅 变化 的 广度 来 看 ,第 12 d 扫 描 时 ,振幅 跨度 最 


表 1 土壤 体积 含水 量 


Tab.1 Soil volumetric water content 


探测 时 间 /年 -月 -日 XA FY FB SIRE (m * ns") 土壤 体积 含水 量 /% 平均 土壤 体 含水 量 /% 
0~10 cm 10~20 cm . 20-30 cm 0~10 cm 10~20 cm 20~30 cm 0~30 cm 
2023-02-24 0.036 0.058 0.081 42.296 35.100 34.481 37.3 
2023-02-27 F 0.038 0.061 0.082 27.872 30.091 29.748 29.2 
2023-03-02 Zu 0.041 0.066 0.085 23.617 27.288 24.969 25.3 
2023-03-05 晴 0.045 0.068 0.087 21.770 24.798 20.818 22.5 
2023-03-08 Zu 0.054 0.080 0.089 12.939 14.942 18.813 15.6 
扫描 天 数 一 第 1d 一 第 3d 一 第 6d 一 第 9d 一 第 12d 
(a) 10 cm-EJZ Ath (b) 20 cmt EAk (c) 30 cm 士 层 处 
M, 
P 0.08 A 0.08 去 0.08 
: 0.06 —Ó À 0.06 : 0.06 
go cures 2 0.04 E] 0.04 
B o0 R 0.02 7$ 0.02 
0 0.2 0.4 0.6 0.8 1.0 0 0.2 04 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 
水 平 距离 /m 水 平 距离 /m 水 平 距离 /m 


图 4 雷达 波 速 变化 


Fig. 4 Radar wave velocity variation diagram 
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扫描 天 数 — Fld — 第 3d 一 第 6d — 第 9d 一 第 12d 
(a) 归 一 化 振幅 变化 (b) 时 间 振 幅 变 化 
000 
2000 
1000 
m 
| 器 1900 
n 
—2000 
—3000 
—4000 
0 500 1000 1500 
频率 /Hz 时 间 /ns 
图 5 雷达 振幅 变化 
Fig. 5 Radar amplitude variation diagram 
A2 样 根 信息 及 GPR 根 点 识别 率 
Tab.2 Sample root information and GPR root point recognition rate 
每 次 扫描 GPR 识 别 的 根 数量 / 根 每 次 扫描 GPR 根 的 识别 率 /% 
TE 土 层 深度 FHIR ”实际 埋 根 " 
级 =] H- 
样 根 编号 n Blom 数量 / 根 扫描 天 数 /d 扫描 天 数 /d 
1 3 6 9 12 1 3 6 9 12 
1-8 10 2.0 8 1 0 2 3 4 12:5 0 25 38 50 
9~16 20 2,9 8 5 T 5 8 8 62.5 87.5 62.5 100 100 
17~19 30 4.5 3 1 0 2 2 3 33 0 67 67 100 
合计 19 0 7 9 13 15 32 37 47 68 79 


大 ,为 0 Hz~1500 Hz; 从 振幅 变化 的 幅度 来 看 ,5 次 
扫描 的 归 一 化 振幅 均 在 0~1 之 间 。 如 图 5b, 不 同 扫 
描 天 数 时 雷达 振幅 差异 明显 (P<0.05), 随 探测 时 程 
的 增 大 ,雷达 振幅 也 逐渐 趋 于 平缓 。 探 测 时 程 的 增 
大 即 表示 土 层 深 度 的 增 大 ,由 于 土 层 逐 渐 增 大 , 平 
均 土 壤 含水 量 趋同 ,对 雷达 波谱 的 影响 也 逐渐 减 
小 ,变化 波形 相似 。 
2.3 土壤 含水 量 对 根 点 识别 的 影响 

由 表 2 可 知 ,5 次 雷达 探测 的 识别 效果 有 较 大 的 
差异 。 总 体 而 言 , 探 地 雷达 根 点 识别 率 随 探测 间隔 
天 数 的 增 大 而 减 小 。 

5 次 雷达 探测 扫描 中 , 根 点 识别 率 在 第 12 dH 
描 时 最 高 ,识别 率 高 达 79%; 第 1d 扫 描 时 , 根 点 识 
别 率 为 32%, 是 5 次 扫描 的 最 低 值 。 从 样 根 埋藏 的 
深度 来 看 ,20 cm 深度 处 根 点 识别 效果 最 好 ,5 次 雷 
达 探 测 扫描 时 , 探 地 雷达 根 点 识别 率 分 别 可 以 达到 
62.5% ,87.596 ,62.596 , 100961 100% ; 10 em 土 层 处 ， 
探 地 雷达 根 点 识别 率 最 低 为 0, 最 大 为 50% ; 30 cm 
土 层 处 , 探 地 雷达 根 点 识别 率 最 低 为 0, 最 大 为 


BZ 
b 


100%. 

如 图 6 所 示 , 计 算出 5 次 雷达 探测 扫描 的 根 点 
距离 均 方 根 误差 ,分别 为 26.6 cm、25.6 cm、24.8 cm, 
21.4 cm 和 18.2 cm。 可 以 看 出 探 地 雷达 根 点 距离 均 
方 根 误差 随 扫 描 天 数 的 增加 而 减 小 , 随 平均 土壤 含 


水 量 的 减 小 而 减 小 。 
30 RMSE 40 
— 平均 土壤 含水 量 
25 人 | 35 , 
E 
a 
: 15 as : 
M 20 * 
5 
15 
0 
1 3 6 9 12 
扫描 天 数 /d 


图 6 不 同 扫 描 天 数 的 探 地 雷达 根 点 距离 均 方 根 误差 
Fig.6 Root-point distance error of ground penetrating 


radar with different scanning days 
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3 讨论 


3.1 土壤 含水 量 对 探 地 雷达 探测 波 速 的 影响 

如 图 7 所 示 ,对 5 次 雷达 探测 时 的 平均 土壤 含 
水 量 条 件 下 10 cm、20 cm 和 30 cm 土 层 深度 的 雷达 
波 速 信息 进行 模型 拟 合 。 探 地 雷达 波 速 与 平均 土 
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数 变 化 明显 。 当 土壤 水 分 较 高 且 处 于 较 饱 和 状态 
时 ,反射 系数 较 低 ;相反 ,干燥 的 土壤 反射 系数 则 呈 
现 较 强 的 响应 ””。 当 平均 土壤 售 水 量 为 37.3% 时 ， 
根 点 附近 的 电磁 波 反 射 系数 平均 值 为 0.145, 而 当 平 
均 土 壤 仿 水量 为 15.6% 时 , 根 点 附近 的 电磁 波 反射 
系数 平均 值 为 0.210。 因 此 ,雷达 波谱 振幅 随 着 平均 


二 含水 量 之 间 呈 显著 的 久 相 关 (P<0.05)。3 个 土 层 
深度 上 线性 拟 合 的 尼 分 别 为 0.916、0.988 和 0.91。 
在 3 个 土 层 深 度 上 , 随 着 平均 土壤 含水 量 的 增 大 , 探 
地 雷达 波 速 在 明显 减 小 (图 4) ,造成 这 种 变化 的 原 
因 可 能 是 由 于 第 1d 扫 描 时 , 样 坑 表层 (0~10 em) + 
壤 含 水 量 高 达 42.296% , ifii F ES (20-30 em) EES 
水 量 相 对 较 低 ,为 34.481% ,由 此 形成 的 水 分 梯度 差 
异 导 致 了 这 种 变化 。 丰 寅 的 水 分 会 导致 士 壤 水 分 
的 过 滤 和 再 分 配 的 非 平 稳 过 程 ,从 而 导致 土壤 的 介 
电 常 数 有 很 大 的 变化 下。 在 本 研究 中 ,雷达 天 线 电 
磁 波 传 播 中 遇 到 的 主要 介质 为 土壤 和 根系 ,在 不 同 
平均 土壤 含水 量 条 件 下 ,雷达 波 速 不 同 , 即 是 电磁 
波 所 过 介质 的 介 电 常数 发 生 了 变化 。 周 赎 山 等 ” 
研究 发 现 不 同 频 率 下 计算 的 雷达 波 速 相同 , 即 雷 达 
波 速 仅 与 介质 的 相对 介 电 常数 有 关 , 与 雷达 波 频 率 
无 关 。 这 为 本 研究 在 使 用 单一 的 500 MHz 频率 天 线 
的 实验 条 件 下 ,得 出 平均 土壤 含水 量 与 雷达 波 速 之 
间 呈 负 相 关 的 结论 提供 了 理论 优 证 。 
3.2 土壤 含水 量 对 探 地 雷达 探测 波谱 振幅 的 影响 

探 地 雷达 的 波谱 振幅 受 平均 土壤 含水 量 的 影 
响 。 较 高 的 土壤 含水 量 会 导致 电磁 波 信 号 衰减 加 
速 , 从 而 使 雷达 波谱 振幅 的 频 度 和 幅度 都 趋 于 减少 
和 平缓 ”。 为 了 更 好 地 人 研究 平均 土壤 含水 量 、 根 系 
和 波谱 振幅 之 间 的 关系 ,本 研究 使 用 电磁 波 反射 系 
数 来 拟 合 各 因素 之 间 的 相关 性 。 

如 图 8 所 示 ,在 不 同 平均 土壤 含水 量 条 件 下 ,可 
观察 到 探 地 雷达 识别 的 根 点 附近 的 电磁 波 反 射 系 


“ (a) 10 cm 土 层 


m 0.08 


T 0.05 。 y=—0.0011x+0.0710 F 
É R?=0.916 a 
= 0.04 c 007 
E. 2 
B i$ 0.06 


0.03 


15 20 25 30 35 
土壤 含水 量 /% 


(b) 20 cm 士 层 


15 20 25 30 35 
土壤 含水 量 /% 


土壤 含水 量 的 增加 而 减 小 。 通 过 分 析 认 为 ,平均 土 
壤 含 水 量 的 增加 导致 了 根系 和 土壤 之 间 的 介 电 篆 
数 差 值 增 大 ,从 而 使 根 点 附近 的 电磁 波 能 量 衰减 ， 
反射 系数 变 小 ,而 波谱 振幅 则 呈现 平缓 的 变化 。 电 
磁 波 的 反射 系数 取决 于 界面 两 边 介 质 的 相对 介 电 
常数 差异 , 当 目 标 物 与 周围 介质 的 电磁 性 参数 有 较 
大 差别 时 ,反射 系数 驶 会 较 大 , 即 反射 波 的 振幅 和 
AE TEIE 。Ryazantsev 等 下 研究 认为 ,土壤 水 分 
显著 影响 雷达 电磁 波 的 传播 速度 和 振幅 ,而 探 地 雷 
达 根 点 附近 的 电磁 波 反射 系数 随 着 平均 土壤 含水 
量 的 增加 而 减 小 。 

5 种 平均 土壤 含水 量 条件 下 根 径 与 控 地 雷达 识 
别 的 根 点 附近 反射 系数 的 拟 合 关系 如 图 9 所 示 。 由 
图 9e 可 以 看 出 ,在 平均 土壤 会 水 量 为 37.3% 时 , 根 径 
越 大 ,电磁 波 反 射 系 数 越 大 。 从 图 9a~9d 可 知 , 根 径 
越 大 ,反射 系数 越 小 。 造 成 这 一 现象 的 原因 ,可 能 
是 由 于 实验 预 埋 的 样 根 在 10 cm 和 20 em 土 层 人 处 较 
细 (d=2 cm 和 d=2.5 cm) ,而 在 30 cm 土 层 处 较 粗 (d= 
4.5 cm) , 样 坑 在 静 置 过 程 中 下 层 土 壤 水 分 较 表 层 土 
壤 水 分 燕 散 更 慢 , 仿 水量 更 加 稳定 ,对 雷达 电磁 信 
号 的 衰减 作用 更 加 明显 ,从 而 导致 下 层 直 径 较 大 的 
样 根 反射 系数 变 小 。 因 此 ,在 探 地 雷达 根系 探测 过 
程 中 ,应 该 尽量 保证 不 同 土 层 深度 的 平均 土壤 含水 
量 相对 稳定 且 水 分 差异 较 小 ,以 此 来 提高 探 地 雷达 
的 根系 识别 率 。 研 究 认为 ,在 使 用 探 地 雷达 对 根系 
进行 探测 时 ,使 用 电磁 波 反射 系数 .振幅 值 以 及 其 
他 物理 指标 去 模拟 根 径 .平均 土壤 含水 量 以 及 其 他 
制约 因子 之 间 的 相关 性 是 可 行 的 。 


(c) 30 cm 士 层 
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图 7 不 同 平均 土壤 含水 量 下 雷达 波 速 拟 合 关系 


Fig.7 Fitting relationship of radar wave velocity under different average soil moisture contents 
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= 020 - 等 J 研究 认为 ,降水 较 多 的 季节 中 , 探 地 雷达 能 量 
$ ol8 | 衰减 更 快 ,根系 探测 精度 会 受到 影响 。De 等 加 和 
0.16 + Sun 等 研究 发 现 , 探 地 雷达 根系 探测 中 ,高 土壤 含 
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图 9 不 同 平均 土壤 含水 量 下 根 径 与 探 地 雷达 根 点 反射 系数 拟 合 关系 


Fig.9 Fitting relationship between root diameter and GPR root reflection coefficient under different average soil moisture contents 
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提高 。 本 人 研究 认为 , 随 着 平均 土壤 含水 量 的 增加 ， 
探 地 雷达 的 根 点 识别 率 减 低 。 在 天 线 频率 不 变 的 
情况 下 (500 MHz) ,进行 了 5 次 雷达 探测 以 识别 根 
系 ,结果 发 现存 在 明显 差异 。 由 于 实验 中 样 根 的 含 
水 量 基本 稳定 ,而 5 次 雷达 探测 时 样 坑 0~30 em E 
层 的 平均 土壤 含水 量 依次 为 37.3% .29.2% 、25.3%、 
22.5% 和 15.6%。 因 此 可 以 推 朵 ,平均 土壤 含水 量 的 
变化 导致 了 根系 和 土壤 的 相对 介 电 常数 差 值 发 生 
变化 ,进而 影响 根 点 识别 的 准确 度 。 仅 从 平均 土壤 
含水 量 角 度 考 虑 根 的 识别 率 , 探 地 雷达 在 进行 根系 
探测 时 , 砂 土 优 于 壤土 ,壤土 优 于 和 莫 土 "。 

在 5 次 雷达 探测 扫描 中 ,不 同 土 层 深度 的 根 点 
识别 率 如 图 10a 所 示 。 在 3 种 土 层 深度 上 ,20 om 
层 深度 的 识别 率 最 佳 , 均 超过 60%。Tanikawa 等 
研究 发 现 , 根 系 的 埋藏 深度 对 探 地 雷达 的 准确 性 有 
影响 。 歼 荤 等 ”也 提出 了 类 似 观点 , 即 根系 的 埋藏 
深度 和 水 平 间 距 都 会 影响 探 地 雷达 的 准确 性 。 在 
10 cm 土 层 处 的 雷达 识别 率 相对 较 低 , 可 能 是 因为 
浅 层 地 面 直达 波 的 信号 较 强 ,容易 掩盖 一 些 浅 层 根 
的 有 效 信号 ,从 而 导致 根 点 识别 率 较 低 。 同 时 , 频 
率 为 500 MHz 的 天 线 在 浅 层 根系 探测 中 的 分 辩 率 相 
对 较 低 ,这 也 容易 导致 漏 探 测 到 一 些 浅 层 根系 , 进 

影响 浅 层 根 的 识别 率 。 

在 30 cm 土 层 处 ,第 1d 扫 描 和 第 12 d 扫 描 的 根 
点 识别 率 分 别 为 33% 和 100%。 两 次 雷达 扫描 时 ， 
0-30 cm 土 层 的 平均 土壤 含水 量 分 别 为 37.3% 和 
Ce 
黏 粒 对 雷达 波 的 干扰 较 强 ,会 加 速 雷达 电磁 信号 的 
SOR. MTT RU, FET BLE UD 
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分 的 蒸发 和 下 渗 , 低 平均 土壤 含水 量 下 雷达 信和 号 的 
衰减 变 缓 ,识别 率 逐 渐 提 高 。 此 外 ,在 第 3 d 扫 描 
时 , 探 地 雷达 在 3 个 土 层 深度 上 的 识别 率 分 别 为 
0%、87.5% 和 0%。 

实验 发 现 , 两 次 扫描 间隔 期 间 , 由 于 水 分 下 渗 
造成 了 样 坑 表层 土壤 下 隐 ,在 探 地 雷达 探测 的 测 线 
路 径 下 方 存在 空洞 ,扫描 后 我 们 对 表层 土壤 下 陷 区 
进行 了 回填 ,以 保证 后 续 实 验 的 顺利 进行 。 因 此 ， 
人 研究 认为 这 一 情况 导致 了 10 cm 土 层 处 雷达 反射 波 
言 号 发 生 剧 烈 变化 ,从 而 掩盖 了 该 土 层 处 的 样 根 反 
射 信号 ,在 后 期 雷达 影像 解 译 过 程 中 无 法 判定 根 点 
人 位置。 同时 ,第 3 d 扫 描 时 ,30 cm 土 层 处 的 土壤 含 
ds M MP E id 

雷达 信号 衰减 的 影响 ,可 能 高 于 其 他 土 层 深 度 。 
ee EL 
以 及 土壤 含水 量变 化 是 造成 第 3 d 扫 描 时 ,10 cm 和 
30 cm 土 层 处 的 根 点 识别 率 为 0 的 主要 原因 。 建 议 
在 进行 相关 实验 时 ,夯实 探 地 需 达 测 线路 径 下 的 土 
EE ,避免 土壤 下 陷 影 响 实验 精 度 。 

在 不 同 平均 土壤 含水 量 条 件 下 , 根 径 变化 对 探 
地 雷达 根 点 识别 的 影响 见 图 10 b。 平均 土 壤 仿 水量 
随时 间或 不 同 地 点 的 变化 会 影响 信号 的 衰减 和 传 
播 速 度 , 从 而 影响 探 地 雷达 返回 信号 的 强度 WY。 在 
2 cm 和 2.5 cm 的 平均 根 径 级 别 中 , 根 径 越 大 , 探 地 雷 
达 的 根 点 识别 率 就 越 高 。 但 是 在 平均 根 径 为 4.5 cm 
时 ,5 个 平均 土壤 含水 量 梯 度 下 的 识别 率 都 小 于 根 
径 为 2.5 em 时 的 识别 率 。 这 可 能 是 由 于 短 时 间 内 灌 
入 样 坑 的 水 分 ,在 30 cm 土 层 处 保持 较 高 且 较 为 稳 
定 的 状态 ,过 量 的 水 分 含量 导致 识别 率 降 低 。 如 果 
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图 10 不 同 影响 因素 下 探 地 雷达 根 点 识别 结 


Fig. 10 GPR root point identification results under different influencing factors 
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考虑 雷达 天 线 的 影响 , 浅 层 根 系 配 合 高 频 天 线 在 低 
平均 土壤 含水 量 和 粗 根 条 件 下 容易 被 探测 到 ,深层 
根系 配合 低频 天 线 在 低 的 平均 土壤 含水 量 和 粗 根 
条 件 下 容易 被 探测 到 1。 
4 结论 

本 研究 通过 控制 实验 ,选用 500 MHz 探 地 雷达 
在 山西 吉 县 获 家 川 流 域 的 油 松 人 工 林 内 开展 野外 
预 埋 根系 的 探测 实验 ,在 5 个 土壤 含水 量 梯度 下 对 
探 地 雷达 波 速 .振幅 等 均 进行 了 分 析 ,结合 探 地 雷 
达 波 谱 图 像 , 证 实 了 土壤 含水 量 对 探 地 雷达 根系 探 
测 有 显著 影响 。 研 究 表明 , 探 地 雷达 波 速 .振幅 、 反 
射 系 数 、 根 点 识别 率 都 与 土壤 含水 量 呈 负 相 关 , 随 
着 土壤 含水 量 的 增加 , 探 地 雷达 波 速 会 逐渐 减 小 ， 
振幅 变化 趋 于 平缓 , 根 点 识别 率 降低 。 在 良好 的 探 
测 条 件 下 , 探 地 雷达 能 够 较 好 的 识别 地 下 粗 根 。 土 
壤 含 水 量 为 13.6% 时 , 根 点 识别 率 最 高 (79% ) ,土壤 
含水 量 为 37.3% 时 ,识别 率 最 低 (32%) 。 因 此 ,研究 
认为 ,土壤 含水 量 是 影响 探 地 雷达 根系 探测 精度 的 
主要 因素 ,土壤 含水 量 较 高 的 土 层 , 根 点 识别 率 相 
对 较 低 。 综 上 ,建议 在 土壤 含水 量 较 低 的 地 区 或 季 
节 ,使 用 探 地 雷达 进行 植物 根系 探测 ,同时 可 以 采 
用 多 种 天 线 频率 组 合 进行 探测 ,从 而 获得 更 高 的 精 
准 度 和 更 有 参考 价值 的 植物 根系 生物 量 研究 结果 。 
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Effect of soil moisture content on the accuracy of root configuration 
detection by ground penetrating radar 


WANG Bo', ZHANG Jianjun", LAI Zongrui, ZHAO Jiongchang, HU Yawei, 
YANG Zhou, LI Yang, WEI Zhaoyang 
(1. School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China; 


2. Jixian National Forest Ecosystem Observation and Research Station, Chinese National 


Ecosystem Research Network, Jixian 042200, Shanxi, China) 


Abstract: Root biomass is an important part of the soil ecosystem, however, due to the limitations of measure- 
ment techniques and methods, it is impossible to rapidly assess it. The ground penetrating radar (GPR) is an effi- 
cient and nondestructive geophysical tool through which root information can be obtained without damaging the 
soil environment. However, accuracy during the detection and identification of underground roots by GPR is sig- 
nificantly affected by many factors such as soil water content, root roughness, length, and extension direction. In 
particular, soil water content has an obvious effect on root detection. In this study, in situ root embedding detec- 
tion experiments were carried out to investigate the influence of soil water content on root detection via GPR. 
Combined with the changes in wave velocity, amplitude, and the root reflection coefficient of the GPR, the root 
point identification rate and root point distance root point root mean square error were analyzed under different 
average soil water content. The results showed that (1) the wave velocity and amplitude of the GPR were impor- 
tant parameters to determine variations in soil water conten during root detection; (2) the GPR' s velocity de- 
creased and the radar amplitude flattened with the increasing soil water content; (3) as the root diameter increased 
at different soil depths, the GPR' s wave velocity also increased and the GPR' s amplitude tended to be drastic 
changes; (4) the root point recognition rate and soil water content were negatively correlated (P<0.05), and the 
best recognition effect was achieved when the soil water content was 15%-25%. These results are great signifi- 
cance for quantifying the effect of soil moisture content on the root detection accuracy of the GPR. In addition, it 
provide a reference for using this tool for root detection and are particularly important for the estimation of plant 
root biomass. 
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